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' In this article we develop Poincare ideas about a heat balance of ideal gas considered as a collisionless continuous medium. 

We obtain the theorems on diffusion in nondegenerate completely integrable systems. As a corollary we show that for 
any initial distribution the gas will be eventually irreversibly and uniformly distributed over all volume, although every 
particle during this process approaches arbitrarily close to the initial position indefinitely many times. However, such 
individual returnability is not uniform, which results in diffusion in a reversible and conservative system. Balancing of 
pressure and internal energy of ideal gas is proved, the formulas for limit values of these quantities are given and the 
| classical law for ideal gas in a heat balance is deduced. It is shown that the increase of entropy of gas under the adiabatic 

■ extension follows from the law of motion of a collisionless continuous medium. 

d 



1. Heat balance 

J> \ The establishment of heat balance of gas in a vessel is one of the central problems of nonequilibrium 
(N 



statistical mechanics. The conventional model is Boltzmann -Gibbs gas: ensemble of a large (but 
finite) number of identical solid balls elastically colliding with each other and with the walls of the 
vessel. According to the classical approach based on the Boltzmann kinetic equation the process starts 
with practically instantaneous establishment of the Maxwell velocity distribution, and then (not so 
<^ • fast and with oscillations) the gas density becomes balanced [Q. 

Unfortunately, such approach involve several fundamental difficulties. First, the Boltzmann equa- 
tion is approximate. It does not take into account multiple collisions, and, besides this approach as- 
sumes the statistical independence of the number of pair collisions. This assumption (Stosszahlansatz 
. by P. Ehrenfest and T. Ehrenfest [2]) is plausible, but it certainly does not follow immediately from the 
dynamics of the Boltzmann -Gibbs gas model. Furthermore, there are difficulties in an adjustment 
$_i ■ of the solutions of the Boltzmann equation with the reversibility property of the dynamics equations 



and with the Poincare theorem on returning (see [31 E] f° r t ne discussion of these problems) . 

Logical opportunity of the adjustment of the irreversible behavior of a system with the properties 
of reversibility and returnability was shown by M. Kac El on a so-called circular model, which, 
however, is not related to the gas theory. A phase space in the Kac model is an ensemble of white and 
black balls in vertexes of regular ra-gon. Besides, a set M of vertexes of n-gon consisting of m < n/2 
elements is selected. The dynamics of the circular model is determined by rotation on one element 
counter-clockwise. If a ball does not belong to the set M, its color does not change, and if it belongs 
to M, the color of the ball changes to the opposite. The dynamics of such system is clearly invariant 
relatively to the direction of rotation (reversibility). It is also possible to verify that after In rotations 
the system will turn into the initial state (returnability). 

Let N c (t) and Nb(t) be the numbers of black and white balls in integer moments of time t. As 
Kac has proved, the mean value of the ratio 

N c (t)-N b (t) x 
n /' 
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calculated over all possible states of the set M, decreases for n — > oo as (1 — 2/x) , where /x is a limit 
value of ratio m/n. Thus, if fi < 1/2 then after the large enough period of time the average number 
of white and the average number of black balls will coincide. 

The Kac circular model is an advanced version of earlier Ehrenfests' model 0, that possessed 
only the property of reversibility. The model also shows a characteristic feature of the conventional of 
reasoning in statistical mechanics: evaluation of average values, passage to the limit on the number of 
particles (n — > oo), and then passage to the limit on time (t — > oo). The latter limit is connected with 
the fact that the average returning time tends to infinity together with n, hence one should choose a 
time interval, which is less than the order of this quantity. 

Actually, the justification of thermodynamics involves additional difficulties of a different kind. 
The matter is that the ideal gas is considered as a system of noninteracting particles. In particular, 
they cannot collide with each other. It is under this assumption that the perfect gas law is deduced 
in statistical mechanics. On the contrary, when the interaction is taken into account (in particular, 
assuming the possibility of collisions is ), then using the canonical Gibbs distribution we obtain the 
equation of state, different from the classical Clapeyron equation (see (0]). 

On the other hand, as it is shown in papers |Hj, the Clapeyron equation is deduced from 
the general principles of statistical mechanics under the assumption that the density of distribution 
of probabilities is a single- valued function of the total energy of system of particles. It should be 
emphasized that this function does not necessarily coincide with the density of Maxwell distribution. 

The ideal gas is the fundamental model of mechanics of a continuous medium and statistical 
mechanics. Therefore, the problem of justification of irreversible behavior of ideal gas that does not 
require the Boltzmann mechanism of pair collisions gains a special importance. And it is not at all 
obvious that such irreversibility mechanism actually exists. 

These problems constitute the subject of the present paper. 

2. The ideal gas as a collisionless continuous medium 

The above mentioned difficulties can be overcome if the ideal gas is considered as a collisionless 
continuous medium. It is necessary to note that the hypothesis on the continuity property of gas is in 
good agreement with the continuity of velocity distribution of gas particles. 

Besides, collisionless models play an essential role in many parts of mathematical physics. As 
an example, the theory by Ya. B. Zeldovich could be mentioned that explains an occurrence of inho- 
mogeneities in a distribution of pulverulent substance in the Universe (see *9J. Another important 
example is the Burgers equation. It describes the dynamics of fluid without pressure and is one of 
possible simplifications of the Navier- Stokes equations [10] . Multidimensional hydrodynamics of in- 
variant manifolds of Hamiltonian systems, developed in work (11 j . also describes the evolution of a 
collisionless medium. 

For the first time, the ideal gas has, apparently, been considered as a collisionless continuous 
medium by H. Poincare in ^2]- He studied the behavior of ideal gas in a rectangular parallelepiped 

IP = {0 < z x < h, < z n < U 

(l s being the edges of the parallelepiped) for n = 1, 2 and 3. Poincare called such gas one- dimensional. 
In Poincare's terminology, three-dimensional gas is consisted of molecules, which can collide with 
each other. His basic observation was that, independently of the initial distribution, gas eventually 
tends to uniform filling of IT. Thus, the ideal gas shows the irreversible behavior. Every particle 
of gas approaches arbitrarily close to the initial position infinitely many times. However, because 
of nonuniformity of the returnability property, a nonreversible diffusion of gas occurs. Besides, the 
equations of motion of a collisionless medium are invariant under reflection of time 1 1— ► —t. Thus, as 
far back as in 1906 Poincare showed on the simplified model (directly related to the kinetic theory) the 
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compatibility of the reversibility and returnability properties with irreversible behavior of a dynamical 
system. 

Unfortunately, these remarkable ideas of Poincare were not properly understood and remained 
unclaimed. I have not found any work on statistical mechanics that mentions his ideas in connection 
with the problem of irreversibility. The comments to the Poincare's work of 1906 in volume III of the 
Russian edition of his collected works completely miss the point. 

Some interesting works on the kinetic theory, that also use the model of collisionless medium, 
have appeared recently. As an example we shall mention the dynamical demon of Maxwell [l'dl I14| . 
However, they do not refer to the Poincare's pioneer works. 

Of course, Poincare has considered only the most simple variants and his works do not contain 
precise statements with complete and rigorous proofs in modern understanding of these words. How- 
ever, his ideas finally (after almost 100 years) deserve an involvement into the area of study of modern 
nonequilibrium statistical mechanics. The purpose of the present paper is the development of the 
Poincare's ideas on heat balance of ideal g collisionless continuous medium. 

So, we consider the dynamics of particles in a ra-dimensional parallelepiped Ii n . Clearly, IF 1 allows 
a natural 2 n -leaf covering by n-dimensional torus T™ = {x\, . . . , x n mod 2tt} with branching on the 
boundary of Tl n . Variables x and z are related as follows: x = ttz/1 if z increases from up to I, and 
x = 2tt — irz/l if z decreases from / to (Fig. 1). 




Fig. 1 

Let v±, . . . , v n be the velocity components of a gas particle in H" C M n = {z}. Then the rates of 
variation of its x-coordinate are equal to 

vrvi irv n , , 

UJ X = —,..., LO n = ——. (2.1) 

Hence, in variables x mod 2tt, u, the dynamics of gas particles in n is described by the equations 

x s = uj s , uj s = (s = l, ...,n). (2.2) 



3. The first theorem on diffusion 

Equations (2.2) describe an evolution of the integrable system. A phase space 

F 2n = p n xJ n = % mQd ^ 

is foliated on invariant tori u = (uj\, . . . , u n ) = const, which are filled with conditionally periodic 
trajectories with frequencies lo\, uo n . For almost all uj £ M. n these trajectories are everywhere 
dense (and, actually, uniformly distributed) on torus u = const. 
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Such picture is generally characteristic for completely integrable Hamiltonian systems with com- 
pact energy surfaces |15j . In a neighborhood of invariant tori it is possible to introduce the action-angle 
variables x mod 2ir, y; in these variables the Hamiltonian equations take the form: 

y s = 0, x s = io s {y); l^s^n. (3.1) 

In a nondegenerate case, when 

d(u>i, ...,uj n ) ^ Q 
9{yi, ■ ■ ■ , y n ) 

it is possible to change from variables y to new variables to. In these coordinates equations (3.1) have 
form (2.2). Thus, equations (2.2) represent a universal form of equations of motion of nondegenerate 
completely integrable systems. 

Let f(uj,x) be a Lebesgue integrable function, which is 27r-periodic on each of coordinates 
x\, x n . Let g: T n — > K be a Riemann integrable function; in particular, it means that it is 
restricted. Let's introduce the following function of time 



K(t) = J J f(co, x-tot)g(x)d n xd n uj. (3.2) 

R n T™ 

Since the function f(u, x — cot) is Lebesgue integrable on all values of t, and g is measurable and 
restricted function, integral (3.2) is correctly defined. 

The function K(t) has an obvious interpretation. First of all we note that it follows from equa- 
tions (2.2) that x — tot = xq = const. Let / ^ be (according to Gibbs) the density of distribution 
of integrable systems in P (density of probability measure), and g be the characteristic function of 
Jordan measurable domain D on T". It is clear that 



(/) = j fd n xd n u; = l. 



It is easy to see, that in this case K{t) is equal to a fraction of all the systems whith phases (i. e. the 
x-coordinates of points on T n ) belonging to the domain D. 
Let's study the behavior of function K(t) for t — > looo. 



Theorem 1. There exists 



lim K(t) = (f)g, (3.3) 

t— >±oo 



where 

T" 

Let's return to the case, where / is the density of distribution of probability measure, and g is 
the characteristic function of measurable domain D. Then relation (3.3) transforms into the equality 

lim Kit) = mesG n . (3.4) 

t^ioo v ' mesT n v ' 

Hence, independently of the initial distribution after an unbounded period of time the system becomes 
evenly distributed on phases. This result shows the irreversible diffusion in nondegenerate integrable 
systems. 

We divide the proof of theoremnonto several items. Since a Lebesgue integrable function can be 
presented as a difference of two non- negative integrable functions, we shall assume / ^ 0. 
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1) Let g(x) = g = const. Then 

J f(u, x - cvt)g d n x d n uj = (f)g 
? 

using the formula of change of variables in a multiple integral. 

2) Let g(x) = expi(m, x), m £ Z n , and By setting u = x — ujt, we obtain 

K(t) = J f^y^d^, 

IE" 

where 

f-m = J /(«>, u)e^d n u 

is a Fourier coefficient of function /, considered as a function on T n , multiplied by (2-7r) n . The Fubini 
theorem implies that the function /_ m : R n — > R is Lebesgue integrable. Hence (according to the 
theory of Fourier transform), K{t) — ► for t — > oo. 

3) Items 1 and 2 imply that theorem^ is valid for any trigonometric polynomial g. 

4) We shall now use the well-known statement from the theory of Riemann integral (compare 
with Let function g: T n — > K be Riemann integrable. Then, for all e > 0, there exist two 
trigonometric polynomials g\ and g<i such that 

(a) 9l (x) ^g(x) ^ 92 (x)forallxGT n , 

(b) g 2 -9i< £■ 

The proof of this statement uses the Weierstrass approximation theorem. 

5) Let 

Kj(t) = f f(u, x - L0t) gj {x)d n xd w ; j = 1, 2. 
? 

Since / ^ 0, then for all t, 

K^t) ^ K(t) ^ K 2 (t). 
According to item 3, for t — > oo the difference K2(t) — K\(t) tends to 

(/) (52 -9i) <</>£■ (3-5) 

Here we use property (b) from item 4. 

Thus, inequality (3.5) holds for all t > 71(e). 

6) According to item 3, for all e > there exists 12(e) such that for all t > 12(e) we have 

\Kj(t) - (f)gj\ < e. (3.6) 
On the other hand, using property (a) and inequality / ^ Owe obtain 

(f)9i^(f)9<(f)92- (3-7) 
Inequalities (3.5)-(3.7) imply that for t > max(Ti(e), 12(e)) the inequality is fulfilled 

\K(t) - (f)g\ < 3{f)e B 

These arguments are similar to the proof of the Weyl uniform distribution theorem As a 

matter of fact, Poincare does not give the precise formulation of the statement on the limit uniform 
phase distribution. He considers only item 2 in the special case of n = 1, assuming that the function / 
is continuously differentiable with respect to to. 
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4. Equalization of density 

Let f(uj, x) ^ be the density of distribution of gas particles in a phase space for t = 0, and g(z) be 
the characteristic function of a measurable domain G in IT™. By definition of density, (/) = 1. On the 
other hand, using the explicit formulas for 2 n deaf covering T n — > IP, we obtain 



— 1 / i \ jn mes G ( a -i \ 



Hence, by theorem^ 



lim K(t)- mesG 



*— >±oo mes II ' 

Thus, after some period of time t a fraction of gas particles, which are situated in the do- 
main Gcll, is proportional to the volume of G. So, we come to the following result (formulated 
by Poincare d]): independently of the initial distribution of /, for t — > +oo and t — > — oo the density 
of gas in the vessel II will be irreversibly equalized. Let's emphasize once again that the diffusion of 
a collisionless medium is determined by nonuniformity of returnability of its particles to the initial 
positions. 

Since the difference K(t) — mes G/mes IT tends to zero, fluctuations of density of ideal gas decrease 
unrestrictedly with a course of time. For the Boltzmann-Gibbs gas consisting of finite number of 
particles, it is necessary, at first, to give rigorous definition of density equalization. Unfortunately, 
no precise theoretical results were presented for this subject so far. However, in any case (because of 
the Poincare theorem on returning), equalization of densities of the Boltzmann-Gibbs gas should be 
accompanied by sustained fluctuations. Problems of numerical simulation of Boltzmann-Gibbs gas 
are discussed, for example, in paper [T7| . 

To estimate the rate of equalization of density of ideal gas, we shall consider a simple example. 
We assume that the velocities of gas particles are subjected to the normal distribution: 

^jMr 6, (4 - 2) 

2 2 , I 2 

Here A is the non- negative measurable function on T n , and 

J \{x)d n x = 1. (4.3) 

T" 

If we interpret (4.2) as the density of Maxwell distribution then the variance a 2 is in proportion to 
Kelvin temperature r. 

In the case in question we have the following equality: 

oo _ <£ 

K(t) - SesG = V' 5m A_ m / -JL^-e****^ (4.4) 
mesI1 ™ J T 



where g m is a Fourier coefficient of the lifting of function g(z) on T n , and (2-7r) n A m is a Fourier 
coefficient of function A: T n — > M. By formula (4.3), we obtain a simple estimate |A m | ^ 1. On the 
other hand, using formula (4.1) we have the following inequality: 



I < mes G 
mes 11 
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Hence, from (4.4) we obtain the inequality 

cr 2 t 2 (m,m) 



K{t) _mes± ^mes±y e ~ 2 (45) 
mesll mesll ^ 

The series on the right hand side converges for all t > (if <r/0) and its sum tends to zero extremely 
fast, when t — > ±00. 

The sum of a majorizing series 

a 2 t 2 (m,m) 

e " ~ 2 ( 4 - 6 ) 

is expressed in theta-functions. It is equal to 

[O3, (0, q)f - 1, 

where q = exp(— o~ 2 t 2 /2), and the third theta-function is defined by the series 



00 



-00 



It is clear, that 9s(v, q) — > 1 for q — > 0. 

Estimate (4.5) is universal in a sense that it does not contain the function A. In particular, it is 
possible to use Dirac's 5-function as A. In that case, the gas at the initial moment is concentrated 
in one point (by the way, this situation does not contradict the hypothesis about collisionlessness of 
the medium). Since a 2 is proportional to r, series (4.6) actually depends on the combination tyfr. 
Therefore, the duration of density equalization process decreases with the growth of temperature 
as 1/yfr. 

Note that in contrast to Kac's model (and to the conventional knowledge about the mechanism of 
heat balance in gases) the equalization of density occurs without preliminary averaging with respect 
to the states and without determination. 



5. The second theorem on diffusion 

Let functions /, g : M n xT" — > M. are integrable together with their squares (with belong to class L2 { 
It is clear that for all values of t the function f(co, x — cot) also belongs to L 2 . Therefore, the function 

K{t) = j f(uj,x-ujt)g(uj, x)d n xd n oj. 

v 

is correctly defined. 

Theorem 2. Under the above assumptions we have 



lim K{t) = (27rf [ fgd n u;. (5.1) 
^±00 J 



This result, certainly, does not follow from theorem^ (as well as theorem ^ is not a corollary of 
theorem [2|. In Poincare's paper |12j formula (5.1) is not mentioned. 

Before we prove theorem we shall make one auxiliary statement. Let 



E 



fMe^ and ^^M^ ( 5 - 2 ) 
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be Fourier series of functions / and g for a fixed value of uj. The Fubini theorem imply that these 
series are defined for almost all u> € M n . Moreover, for almost all uj the functions / and g belong 
to L 2 {T n ). Hence, 

J f(uj, x - ujt)g(uj, x) d n x = (2vr)" £ /^-^e"*^*. (5.3) 

Since /, g £ L 2 , then the functions \f m 9-m\ are integrable in R n for all m £ Z n . Let's denote 
g'm = 5mexp[z(m, w)*]. It is clear that g' m g'_ m = g m g- m . 



Lemma 1. 



JUmS'- m + fmS'J<T l ^< /(/ 2 +9 2 )<i"lrf'V. (5.4) 



It is a variant of the Bessel inequality. The inequality shows, in particular, that the series in the 
left-hand part of (5.4) converges uniformly on t £ M. 

Proof of the lemma. We shall use an obvious inequality 

\frri9—m f—m9m\ ^ fmf—m 9m9—m- 

Since f m and /_ m are complex conjugate, fmf—m 0. Similarly, 9m9—m 

> 0. 

Now we are to prove the inequality 

J2 j fm.f-md n UJ^ f fd n xd n 0J. (5.5) 

m R n P 

The similar inequality holds for function g. 

Indeed, let /jv be the finite sum of the terms of Fourier series (5.2), such that \m\ < N. Then 



0^ / "(/ - f N fd n xd n uj = J f 2 d n xd n oj- j 



fmf—m d LO. 



ml <A r pn 



Hence, inequality (5.5) is valid for any finite sum of the series in the left-hand part of (5.3). We obtain 
the required statement by passage to the limit for \m\ — ► oo. 

Now, let's prove theoremEl By lemma^ (using the well-known theorems by Levy and Lebesgue) 
series (5.3) converges for almost all oj, and it is possible to integrate the series term-by-term. Inte- 
grating both parts of equality (5.3) over W 1 , we obtain the relation 

K(t) = (2n) n [ Jgd n oj + (2vr)" V / f m (uj)g„ m ^)e- lMt d n uj. 



Since functions f m 9-m are Lebesgue integrable, each term of series in the right-hand side tends 
to zero, when t — > ±oo. According to the lemma, this series converges uniformly on t. Therefore for 
any e > there exists N(e), such that the sum of terms of series with the indices \m\ > N(e) is less 
than e/2 for all values of t. The finite sum of remaining terms tends to zero when t — > ±oo. Hence, 
there exists T(e) (actually depending on N(e)), such that for \t\ > T(e) this sum will be less than e/2. 
So, for \t\ > T{e) the sum of series is less than e, q.e.d. 

Theorem |21 make it possible to solve the problem of evolution of density of distribution of / for 
t — > oo. At first, the density f(u, x — ut) seems to oscillate conditionally-periodically and, therefore, 
there is no limit for t — > oo. However, the density of distribution of probabilities does not "exist" by 
itself, but only as an averaging of some fixed function from L 2 . Therefore, the evolution of / for t — > oo 
should be considered in the generalized sense, as it is usually done in the theory of the generalized 
functions (see, for example, [T5|). 
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To understand, how integrable system (2.2) is distributed in phase space P for t — > oo, we shall 
introduce the characteristic function g of the following set 

G = {x mod 2ir, uj: x s ^ x s ^ x s " , u' s ^ lo s ^ u; s ", 1 ^ s ^ n}. 

It is clear that 

( n 
(2vr)- n Y[(x' s - x s ), if J ^ w < w", 



(.9) 



for the remaining lj. 



It is asserted that for t —> oo function f(u>, x — Lot) converges weakly to /. Limit density / is an 
integrable non-negative function on R n , and (/) = 1. 
Indeed, by theorem |2 for t — > oo 



/(w, x - wt)g(u, x) d n x (Tlo -» (2vr) n J fgd n uj 
f\(x" -x' s ) j ' ... J fdhJx... dw n . 



s=l 

But exactly the same result is obtained by direct computation of the average value of limiting density / 
over domain G. 

As an example, we consider a nondegenerate Hamiltonian system with one degree of freedom 
and show that any function of distribution from Li converges weakly to the function depending only 
on the total energy. Let's recall the definition of nondegeneracy. We assume that the whole phase 
space consists of finite number of pieces, invariant relatively to the phase flow; on each of such pieces 
it is possible to introduce action-angle variables y, x mod 2n. The transition from usual canonical 
variables p, q to variables x, y is a symplectic transformation: its Jacobian is equal to one. In new 
variables, Hamiltonian H(p, q) depends only on y. We call the system nondegenerate if 

dy 2 

on each of the invariant pieces. It is possible to verify, for example, that a usual pendulum in the 
gravity field satisfies these requirements. 

In the action-angle variables the Hamiltonian equations 

. = _dH_ ■ = dH_ 

dq ' dp 

become 

x = uj{y), y = 0, (5.7) 
where u = -j— is the frequency of periodic motion. 

According to (5.6), ^ ^ 0. Hence, w is a monotonic function of y and, consequently, it is 

possible to change from variable y to frequency oj on each invariant piece. Then, equation (5.7) takes 
the universal form 

x = u, Cj = 

and we can use theorem |2J the initial density / tends to the "equilibrium" limit /, depending only 
on u>. By frequency, the limit density is a function of y and, hence, of H. 
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The attempts to prove that for t — > oo the distribution of probabilities tends (in some sense) to 
the stationary state that corresponds to the heat balance ^§1 (chapter XII) can be traced as far back 
as to Gibbs. According to M. Kac [1] (chapter III), the idea that probability should be introduced into 
mechanics only by means of the initial density, certainly seems to be very attractive. But, generally 
speaking, this point of view is, apparently, unfounded and the probability should also appear in 
mechanics by other ways. 

In our opinion, theorems H an d El on diffusion in integrable systems show the fruitfulness of Gibbs' 
approach and indicate that his ideas have not yet been fully realized. 



6. Pressure, internal energy and the equation of state 



Theorem U establishes the law of equalization of density of ideal gas in a rectangular parallelepiped. 
Theorem [2] let us prove the equalization of pressure and density of energy of ideal gas, specify the 
formulas for the limit values of these quantities and, thus, deduce the equation of state of gas in heat 
balance. 

First, let's deduce the formula for the pressure of gas on one of the walls from collisions of the gas 
particles with this wall. The deduction of the formula for the pressure follows the classical reasoning 
in the elementary kinetic theory of gas (see, for example, [20]; it is usually assumed that the gas is 
already in heat balance and the velocity distribution is the Maxwell distribution). 

For the sake of clarity, consider a wall, determined by equation z\ = l\. Let's choose an infinitely 
small rectangle a with the center in point z 2 , z 3 ), its area being equal to da = dz2dz 3 . The 
particles of gas that can hit rectangle a with velocity lo = ui2, ^3) in the moment of time t were 
situated in the moment of time t — dt on the parallel rectangle a with the same area da (Fig. 2), 
during period dt they sweep volume 

dv = uj\ dt dz2 dz 3 . 
The number of such particles is equal to 

dn = Nf(u, h, z 2 , z 3 ) dv, 

where the constant coefficient N is equal to "the number of particles of gas 
in the vessel" . Let m be the mass of "one" particle. Then mN is the total 
mass of gas. Numbers m and N are of the conventional character, they are 
introduced for the purpose of comparison of the formulas obtained below 
with the known formulas of statistical mechanics. Since (/) = 1, then mN f 
is the density distribution of the mass of gas. 

During time dt the particles deliver the following impulse to the wall 




dP = 2mwi dn = 2mNuJ\f dtda. 



Fig. 2 



It is clear that dP/dt is a force of pressure. If we divide it by the area da, 
then we shall obtain the elementary pressure on the wall in the point with coordinates Zi, Z2, z 3 . 
Integrating over all velocities, we shall obtain the total pressure 



00 00 00 



P 



2mN 




wjf(u>, h, z 2 , z 3 )du}idL02du) 3 . 



(6.1) 



Let / be an even function of velocities: f(—u>, z) = f(u, z). Then formula (6.1) becomes more 
symmetric: 

p = mN J uj\J{uj, li, Z2, z 3 )d 3 uj. (6.2) 
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Let's emphasize, that in the initial moment of time the pressure p depends on a point on a wall 
z\ = l\. The formula for pressure on the other walls is the same as (6.2), with only change of lo\ 
with lo\ and io\ accordingly. 

Remark 1. How shall we determine the pressure of gas in an arbitrary interior point z € II on small 
surface with a normal vector nl If particles of gas hit the surface from the opposite sides, then (in accordance 
with the assumption of evenness of the distribution function) we shall, obviously, obtain zero pressure. If we 
put a body of small volume into the vessel II, then the particles of gas will hit only one side of the surface. 
Hence, at this moment of time the pressure can be determined by the slightly improved formula (6.2): 

p = mN J (w, n) 2 f(uj 7 z) d 3 u>. 

Now let t go to infinity. Then, according to section 5, density / tends to "the equilibrium" state /, 
depending only on w. As a result, the pressure does not depend on the point of the wall any more: 

p = mN I u>ffd 3 w. (6.3) 



Remark 2. Actually formula (6.3) requires additional justification, since in (6.2) there is no averaging 
by configurational space, and function u> \— > certainly, does not belong to class L2. We should understand 
formula (6.3) differently. 

Let f(co, x) be a density function, "lifted" on direct product R 3 xT 3 . Let's assume 



p{t) = mN j u>(f(u>, x - ut) d 3 uj. 

R 3 

Further, let 

k 

be a Fourier series of function /, defined for almost all w £ I 3 . Then 

p = mN J Loffd 3 LU + mNj2 el{k ' X) J ^fk(uj)e- l{k ^ ]t d 3 uj. 

If the functions u/ffk are summable for all k £ Z 3 ,then p(t) tends to (6.3) for t — > 00. 

Note, if / depends on u> 2 = Y2 u s> t nen the Pascal's law is valid: the pressure in all directions is 
identical. In the opposite case it is not so: a fraction of the gas particles moving in different directions 
is not the same. 

For the average kinetic energy of one particle we have the formula 

e = J f ^/(c, z)<$zJ?lo. 
r 3 n 

If we change to c-leaf covering t 3 -» n, we can write this formula in the following form: 

hhh / / muj 2 f ( \ r3 j3 

j[oj, x) a x a u. 



(2vr) 3 J J 2 

E 3 T 3 

The average (interior) energy of the whole gas is, obviously, equal to 



E= mNv j ^Jd'ca, 



»3 
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where v is the volume of the vessel. The integral on the right-hand side converges if the values of 
pressure of gas are finite on each wall of the vessel. 

If / depends on to 2 , then (6.3) and (6.4) imply the simple relation 

E = (6.5) 

This formula is well known in the theory of ideal gas. Let 

Vt = dE + pdv (6.6) 

be a 1-form of heat inflow; it is not a total differential. In accordance with the law of degradation 
of energy of thermodynamics, for f2 there exists an integrating multiplier: j3Vl = dS, where S is 
an entropy, /3 = l/(kr), r is Kelvin temperature, k is the Boltzmann constant. According to (6.5) 
the integrating multiplier for form (6.6) is equal to (pv) -1 . Hence, pv is proportional to Kelvin 
temperature r, and we obtain the perfect gas law (the Clapeyron equation). 
It is possible to obtain this result differently, assuming 

7 = cp(ftz£), (6.7) 
where /3 = l/(kr), k is the Boltzmann constant, c is a normalized multiplier: 



l = vl p(p™f)d?u. (6. 



The distributions of form (6.7) were studied in paper jHJ; the classical Maxwell distribution also belongs 
to this type. 

After a change uj s = uj s /y/m(3 we obtain 



»3 



J=V(f)^ 

R 3 



P NkT 

h = x — - — , x 



For the Maxwell distribution x = 1. In paper [S] a class of non-Maxwell distributions was described, 
for which the equality x = 1 also holds. 

7. Entropy 

As it is known, the entropy is determined by equality 

S = - J f(u, x) In /(w, x) d n x d n oj, (7.1) 
? 

where / is a function of distribution of probabilities. For the Maxwell distribution it coincides with 
the notion of entropy used in the equilibrium thermodynamics. 

The question is, what is the evolution in time of the entropy in the case in question, when the gas 
is represented as a collisionless continuous medium. It is important to emphasize, that the evolution 
of state of gas is an adiabatic process: there is no energy transfer. 
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To express S as a function of time, it is necessary to substitute x with x — cut in the integrand 
of (7.1). However, for such substitution integral (7.1) does not change and S as a function of time is 
constant. 

This simple observation corresponds to the Poincare's result ^2] that the fine entropy of mathe- 
maticians, in contrast to the rough entropy of physicists, is always constant. By the way, the division 
of the entropy into the fine and the rough corresponds essentially to the fine-grain structure and 
coarse-grain structure of the phase space, that was introduced by T. Ehrenfest and P. Ehrenfest in 
their well-known work [2]. 

It is possible to approach the problem of behavior of the entropy from the other side. We already 
saw in section 5, that the function of distribution f(uj, x — cut) for t — > oo tends in the generalized 
sense to the average value /, depending only on uj. Let's assume that 

Soo = -(2tt)™ f JlnJdTu. (7.2) 



This expression can be interpreted as an entropy in steady equilibrium state. By the way, formula (7.2) 
can be obtained by theorem |5J using expression g = In/. 
The following inequality 

S s$ Soo, (7.3) 

expresses the law of degradation of energy of thermodynamics for irreversible processes. 

To prove (7.3) we shall fix the value of uj and then denote f(u>, x) as p(x). Now, we establish the 
inequality 

p In p ^ ~p In p. 
It is in turn equivalent to the discrete inequality 

^pilnpi ^ (^2 Pi) ln ^Pi 
for the positive pi, which is the special case of the Jensen inequality for a convex function p i— ► pin p. 

m 

As Poincare has noted ^2]) the values of entropy can be compared only in the states of steady 
equilibrium. 

Let's consider a simple, but instructive example. Let vessel n be divided by a barrier into two 
parts and the gas be initially concentrated in one of the parts of II, and in heat balance. Its entropy 
we shall denote by S-. Now we remove the barrier. The gas will extend adiabatically, uniformly 
filling (by theorem ^) the whole volume of II. Let S+ be the entropy of gas after the establishment of 
heat balance, that will happen after an infinite time. Accordingly to (7.3), <S_ ^ S+. Moreover, it is 
possible to show that this case involves the following simple relation: 

5+ = 5_ + ln^, (7.4) 

where V-(v+) is the volume of Il_(n). 

Indeed, let /_ be a density of distribution of gas, being in heat balance in vessel II_; this function 
depends only on velocity u. Then, obviously, 5_ = —V-J, where 

J= f /_ln/_d* W . 



After removal of the barrier the equilibrium is broken and the density f(u), z) now depends on the 
point z £ IT: in the initial moment / = /_, if z G II _ and / = 0, if z € II \ II_. Theorem |21 implies 
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that for t — > oo density / tends in generalized sense to the average value 

n 

Finally, 

i" 

= -u_ J - t>_ I f_d n uj In ^ = S- + In JJt. 



The formula (7.4) coincides with the well-known formula of the increase of entropy for the process 
of free expansion of gas into the void. However, we have obtained this formula, not basing on the laws 
of thermodynamics, but using only the law of motion of ideal gas as a collisionless medium. 



8. The change of the vessel shape 

Will our deductions change if we replace the rectangular parallelepiped with a surface of arbitrary 
shape? This problem has the principal value, not only from the point of view of thermodynamics. 

The matter is that trajectories of particles of gas are essentially the light rays. Therefore, the 
problem can be reformulated in terms of geometrical optics. Let's put a light source (probably dis- 
tributed) inside the closed reflecting surface. The question is, will the illumination inside this surface 
be constant or will it depend on the coordinates? Actually the related problem also arises if we con- 
sider the radiation in a closed volume with the beam approach (see, for example, [21] )■ The above 
question is traditionally answered positively. But, in this situation together with the arguments of 
dynamical character the requirement of heat balance is usually used. By the way, Poincare himself 
made contradictory statements on this subject (see (22] )• 

Actually the answer is certainly negative, and it is easy to understand, bearing in mind the 
presence of focal points and caustics. As a simple example, it is possible to consider the vessel having 
a shape of ellipse and put the light source in one of the focuses. It is easy to understand, that the 
resulting illumination will concentrate on the transversal line of ellipse (by the way, it is unstable) . If 
the light source is not located in the focus, the illumination intensity inside the ellipse will be variable. 
The existence of a limit distribution in this case (as well as in any other integrable problem) follows 
from theorem |21 If billiard is not integrable, the existence of a limit density of distribution presents 
an interesting problem. 

In some sense, any billiard can be arbitrarily precisely approximated by an integrable billiard. On 
a plane it is a polygon, with angles comparable to tt (see, for example, |23|). Except for the integral of 
energy, such systems allow the integral in the form of velocity polynomial. For example, in case of the 
rectangle the additional integral has the power 2 (squared projection of velocities on any of the sides 
is preserved). It is easy to understand, that the billiards showed in Fig. 3 share the same property. 
Using the integrability of these systems, with the help of theorems ^ and |21 it is possible to study 
diffusion of ideal gas as a collisionless medium in vessels of the indicated shape and, in particular, to 
study the processes of diffusion and mixture of gases in vessels with barriers. 

So, the indicated mechanism of irreversible diffusion of collisionless medium is not universal. 
However, for the rarefied gases it should be considered together with the mechanism of chaotization 
based on molecular collisions. According to Poincare, after a time, sufficient enough to let every 
particle go the whole length of the vessel several times, but short enough for the collisions not to 
be too numerous, the mode will be established in gas that corresponds to the equilibrium state of a 
collisionless continuous medium. But this equilibrium will not be final, the collisions will try to break 
it, and only after a much larger time interval the gas will at last reach the final heat balance. 
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By the way, "a strange kinetics" with participation of the dynamic demon of Maxwell discovered 
in papers |13U14| . is also based on consideration of a collisionless medium. The numerical calculations 
show that the density of gas is not equalized in vessels, made up from scattering billiards. As we just 
have shown, the similar effect occurs for integrable billiards. The strangeness of kinetics disappear as 
soon as we begin to take into account the interaction of particles of gas. 

In conclusion of the paper we shall make two remarks. 

Let's consider a nondegenerate and quite integrable Hamiltonian system, which is described by 
equations (2.2) in each invariant domain D x T n of the phase space, where D is a domain in M n = {ui}. 
Let's perturb slightly the Hamiltonian function. Though the perturbed system is not integrable any 
more, the majority (in the sense of Lebesgue measure) of invariant tori do not disappear, but become 
only slightly deformed (the Kolmogorov's theorem on persistence of conditionally- periodic motions). 
When the Hamiltonian function is sufficiently smooth the dynamics on the invariant set of large 
measure will be described again by equations (2.2), but now u G D\M, where mesM — > 0, when the 
perturbation disappears In this situation theorems ^ and [3 will be still valid, but the density of 
distribution of probabilities / must be equal to zero on the direct product M x T n (in the gap between 
Kolmogorov tori). Thus, it is possible to speak about the diffusion of perturbed Hamiltonian system 
on the invariant Kolmogorov set. Note that D\M has the structure of a Cantor set; in particular, it 
is nowhere dense in D. 

Note also that the limit density of distribution / can be obtained by averaging of function / over 
trajectories of system (2.2). Let's assume that 

T 

f(u, x) = Jim ^ / f(uj, x - ujt) dt. (8.1) 

T^oo 1 J 


The Weyl theorem implies that this limit exists for all phases x and f(to, x) = f(to) for all nonresonance 
sets of frequencies to = (u>i, . . . , u n ). We know that the frequencies ui, . . . , u n are in resonance, if 
k\cu2 + . . . + k n u; n = for certain integer k s , not all of which are equal to zero. Since the resonance 
sets ijj £ M. n amount to a set of zero Lebesgue measure, then, from the point of view of the measure 
theory, functions / and / are equivalent. The function / is continuous for nonresonance values u 
and, in general, is discontinuous on a set of the resonance tori [21] (as a classical example of Riemann 
function, continuous in irrational and discontinuous in rational points of the real axis). 

Let's consider now a dynamical system, more general than (2.2). It is defined by the autonomous 
system of differential equations 

x = v(x, u), Co = (8-2) 

on a direct product P = A x R m , where A = {x} is a compact ?i-dimentional manifold; v is a smooth 
vector field on A with an invariant measure dfi = X(x, uj) d n x: 

^ axi 
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In contrast to (2.2), the dimensions m and n do not coincide. For example, for m = 1 equations (8.2) 
describe the Hamiltonian systems with compact energy manifolds; the total energy of the system acts 
as variable ui. 

Let f(u), x) and g(uj, x) be the functions, summable with a quadrate, defined on the phase space P. 
Following the construction of (3.2), we consider the function 

K(t) = f f(u, g-\x, u)g(u, x)) d m co dp, (8.3) 

i? 

where g l is a phase flow of system (8.2). For the purposes of justification of the kinetics of a collisionless 
continuous medium in the general case, we shall prove the following statement: if for almost all 10 the 
dynamical system x = v(x, ui) is ergodic on A, then 

\imK(t)= [ JgmesAd m uj, (8.4) 

t— *oo / 



where 



( - ) = rniA/ ( - )d/i ' mesA = /V 



A A 

The relation (8.4) contains the formula (5.1) as a subcase. 

Note that if for almost all u> the dynamical system on A has the mixing property the formula (8.4) 
is, certainly, valid. Indeed, for fixed oj (by definition of mixing) 

f(uj, g~\x, u))g(u, x)dfi -> J f dfi J gdfi (8.5) 

A A A 

for t — > oo (see, for example, Then we should integrate this relation over R m . 

It is important to emphasize, that the system x = lo on torus is not mixing. In this case when 
proving of (8.4), we essentially use an averaging on cj and the property of Fourier transform. 

Note that for ergodic systems without mixing relation (8.5) is valid only for Cesaro convergence, 
when the left-hand part is time averaged (as in (8.1)) (see [2Sj)- However, the time averaging (charac- 
teristic for the ergodic theorems by von Neumann and Birkhoff ) should be replaced with the averaging 
by u in our case. Note, that in Kac's circular model there is also an additional averaging over all the 
possible states of the set M , which is used in description of the dynamics of white and black balls. 

Let A = T n and the field v does not depend on x. Then it is possible to set A = 1 and, consequently, 
mes A = (27r) n . Formula (8.4) is obviously valid in this case if we assume that m-dimensional surface 
to i — ► v(uj) in PJ 1 = {v} is transversal to the resonance planes (k, v) = 0, k £ Z n \ {0}. The problem of 
validity of (8.4) in the general case is still open. 

The paper is prepared with the financial support of RFBR (99-01-01096) and the "Leading sci- 
entific groups" grant (00-15-96146). 
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